Towards absolute calibration of optical tweezers 
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Aiming at absolute force calibration of optical tweezers, following a critical review of proposed the- 
oretical models, we present and test the results of MDSA (Mie-Debye-Spherical Aberration) theory, 
an extension of a previous (MD) model, taking account of spherical aberration at the glass/water 
interface. This first-principles theory is formulated entirely in terms of experimentally accessible 
parameters (none adjustable). Careful experimental tests of the MDSA theory, undertaken at two 
laboratories, with very different setups, are described. A detailed description is given of the pro- 
cedures employed to measure laser beam waist, local beam power at the transparent microspheres 
trapped by the tweezers, microsphere radius and the trap transverse stiffness, as a function of radius 
and height in the (inverted microscope) sample chamber. We find generally very good agreement 
with MDSA theory predictions, for a wide size range, from the Rayleigh domain to large radii, 
including the values most often employed in practice, and at different chamber heights, both with 
objective overfilling and underfilling. The results asymptotically approach geometrical optics in the 
mean over size intervals, as they should, and this already happens for size parameters not much 
larger than unity. MDSA predictions for the trapping threshold, position of stiffness peak, stiffness 
variation with height, multiple equilibrium points and 'hopping' effects among them are verified. Re- 
maining discrepancies are ascribed to focus degradation, possibly arising from objective aberrations 
in the infrared, not yet included in MDSA theory. 

PACS numbers: 87.80.Cc, 07.60.-j,87.80.-y 



I. INTRODUCTION 

Optical tweezers have become invaluable tools for mea- 
suring forces and displacements at the single molecule 
level in cell biology Transparent microspheres are 
employed as handles and force transducers. In the usual 
domain of application, the trapping force is well described 
by Hooke's law, so that force calibration amounts to de- 
termining the trap stiffness. Commonly applied methods 
rely on comparisons with fluid drag forces or on the de- 
tection of thermal fluctuation features [2j]. 

Absolute calibration should be based on a theoretical 
expression for the stiffness in terms of measurable pa- 
rameters. In the most widely employed optical tweezer 
setup, an incident TEMoo laser beam of vacuum wave- 
length Ao, typically around 1/im, is focused by a high nu- 
merical aperture oil immersion objective of an inverted 
microscope through a glass slide onto a transparent mi- 
crosphere immersed in water inside a sample chamber. 
One finds that high trapping efficiency requires a beam 
waist that slightly overfills the objective. The radius a of 
the microsphere usually ranges between 0.1 Ao and 10 Ao- 
Thus, an accurate theoretical model should ideally satisfy 
several requirements: 

(i) It should realistically model the strongly focused 
laser beam produced by the objective. Since most micro- 
scope objectives are corrected for the visible, there may 
be objective aberrations to account for in the infrared. 

(ii) With oil immersion, the index discontinuity be- 



tween the glass slide and the water gives rise to spher- 
ical aberration, which degrades the focus and must be 
included in the beam description. 

(iii) The direct interaction between the focused beam 
and the microsphere spans the full range between 
Rayleigh scattering and the ray optics limit, so that it 
should be described in terms of Mie scattering. 

(iv) At the upper end of the range of microsphere radii, 
the results should approach the ray optics approximation. 
However, as is typical of "semiclassical" approximations 
[3[, this is not to be understood as a pointwise limit, but 
rather as a size average over rapid interference oscilla- 
tions. 

(v) The directly scattered beam undergoes multiple 
reflections between the microsphere and the walls of the 
sample chamber (reverberation). Usually, only the glass 
slide may be close enough to have a possibly significant 
effect. 

A brief critical review of proposed theoretical models, 
as far as possible in chronological order, follows. 

A ray-optics model was formulated by Ashkin [1] from 
preliminary results derived by Roosen [5[ . It did not take 
into account that high numerical aperture objectives are 
designed to satisfy the Abbe sine condition, which was 
later incorporated into the model by Gu et al @. 

A widely employed wave-theoretical model is based on 
so-called GLMT (Generalized Lorentz-Mie Theory), de- 
veloped by Gouesbet et al. [7|. The laser beam is rep- 
resented by an attempted improvement on the parax- 
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ial Gaussian TEM o model, including fifth-order correc- 
tions inpowers of the ratio between wavelength and beam 
waist jsf. However, it has been shown [9] that such an 
approximation does not correctly represent the field near 
the focus of a high NA objective. 

A new approach was taken by Tlusty et al. [10]. They 
proposed to approximate the trapping force by the gra- 
dient of the dipole interaction energy integrated over the 
trapped particle, taken with respect to an unperturbed 
incident Gaussian beam around the focal point. This 
leads to a Gaussian fall-off of the transverse stiffness as 
a function of a for a ^> Ao, a result that is incompat- 
ible with the hyperbolic fall-off implied by dimensional 
arguments [ll[ , and which is reflected in the large incon- 
sistency they find with experimental data in this region. 

Already in Ashkin's work [4] it was remarked that the 
proper wave description of a highly convergent beam is 
not a Gaussian, but rather the electromagnetic gener- 
alization [12] of Debye's classic exact scalar representa- 
tion formulated by Richards and Wolf. This representa- 
tion, with proper accounting for the sine condition, was 
adopted by Maia Neto and Nussenzveig [11] to evalu- 
ate the axial trapping force, with scattering described 
by Mie theory. The force at the geometrical focus shows 
rapid near-sinusoidal oscillations, accounted for by a sim- 
ple interferometer picture. The axial stiffness approaches 
the Rayleigh limit for k±a <C 1 (k± = 27rni/Ao is the 
wavenumber in the sample region, of index m) and shows 
attenuated interference oscillations for k\a ^> 1, with size 
average approaching the ray optics result. 

In a new treatment by Rohrbach and Stelzer [l3j |. an 
angular spectrum (Fourier) representation of the incident 
beam [l4] analogous to the Richards- Wolf representation 
is employed. However, instead of the exact Mie theory, 
they extend the approach of Tlusty et al. by splitting 
the trapping force into a gradient force and a scatter- 
ing force. The gradient force is given by an expression 
similar to that of Tlusty et al. (apart from the modified 
incident beam). The scattering force is obtained in the 
Rayleigh-Gans approximation, in terms of extinction and 
scattering efficiencies. While the result would suffer from 
a similar incompatibility with the ray optics limit, this 
problem does not arise, because use of the Rayleigh-Gans 
approximation already restricts the domain of applicabil- 
ity to the range 



2(n 2 - ni)a ~ A 



(i) 



where ri2 is the refractive index of the sphere. 

The MD theory was extended to transverse trapping 
forces by Mazolli et al. [15]. It was explicitly demon- 
strated that the exact partial- wave result approaches the 
ray optics one for large k\a, in the sense of a size aver- 
age, as it should. Besides the transverse trap stiffness, 
the equilibrium position with an external applied force 
and the maximum transverse force were also evaluated, 
taking due account of the interplay between axial and 
transverse equilibrium. As remarked by Merenda et al. 
[l6[ 5 it is the only wave-theoretical treatment where this 



effect was correctly incorporated in the evaluation. Only 
a few experimental points were available for comparisons, 
and discrepancies were attributed to the effects of spher- 
ical aberration at the water/glass slide interface [effect 
(ii) above], not included in the MD theory. 

A recent contribution by Rohrbach [ItJ does not prop- 
erly refer to optical tweezers, but to a proposed new in- 
strument, the photonic force microscope [l8j], which em- 
ploys an upright microscope with a water-immersion ob- 
jective, thereby avoiding this spherical aberration effect. 
The theory is the same as in Ref. [13j, and therefore 
it is also inconsistent with the ray-optic limit, although 
the author asserts that it delivers reasonable results even 
when condition (pQ) is not satisfied. 

An improved version of MD theory, including the ef- 
fects of spherical aberration at the water/glass slide in- 
terface, denoted as MDSA theory and described in Sec. 
II, was experimentally tested and the results were briefly 
reported by Viana et al. [19[. The present work is a de- 
tailed presentation and discussion of these results. 



II. MDSA THEORY 

The effect on a focused beam of spherical aberra- 
tion produced by refraction at the interface between two 
transparent media on a focused beam, extending the 
Richards- Wolf solution, has been treated by Torok et al. 
[20||. We follow the same procedure to incorporate this 
effect into the MD theory. Since the refractive index n of 
the glass slide is larger than the refractive index m of the 
water in the sample chamber, part of the incident beam 
angular spectrum of plane waves may exceed the critical 
angle. However, we shall neglect possible contributions 
from evanescent waves, expected to be negligible at dis- 
tances from the interface larger than the wavelength, as 
will be assumed. Thus, for each plane wave component 
of the incident beam, the transmission amplitude is 



T{e) 



2 cos 



cos 6 + N cos 0\ 



(2) 



where N — ni/n, 6 is the angle between the wavevector 
component k at the glass slide (with k = 27rn/Ao) and 
the z-axis, and 0\ = arcsin(sin#/A/') is the corresponding 
angle in the sample chamber [211 ]. 

The refraction also (and more importantly) modifies 
the phases of the different plane- wave components of the 
laser beam. This is quantified by the spherical aberration 
function [20[ 



9(z,0) = k 



L 

~N 



cos# + N(L + z) cos 6i 



(3) 



where L is the distance between the interface and the 
paraxial focal plane (Fig. [TJ). 

The position of the center of the sphere with respect to 
the paraxial focus is (p, 0, z) in cylindrical coordinates. 
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FIG. 1: Schematic representation of the dielectric microsphere 
(radius a) in the water solution. The paraxial focus F corre- 
sponds to the intersection of paraxial rays. Non-paraxial rays, 
not represented for clarity, intersect the symmetry z-axis be- 
tween F and the glass- water plane interface, which is at a 
distance L from the paraxial focal plane. The microsphere 
position is measured with respect to F. 



For simplicity, we take circular polarization. In this case, 
the optical potential does not depend on (j). With respect 
to the ideal case considered in Ref. [15[, the multipole 
coefficients of the incident beam are modified by multi- 
plying each plane- wave component by T{6)e l ^^ z ^ : 



dO sin V cos exp(- 



2 • 2 

-7 sin 



6)T{6) (4) 



In the absence of the interface (N = 1), A coincides with 
the filling factor 1 — exp(— 2j 2 sin^o) defined in Refs. pj| 
and pj]. 

Additional transmission losses through the objective, 
if any, must be taken into account to evaluate the ac- 
tual power P at the sample (rather than the power at 
the entrance aperture) Investigation of the microscope 
objective employed in the present work has revealed that 
its transmittance at the infrared laser wavelength is not 
radially uniform [23]. This affects not only the power, 
but also the intensity gradients in the sample region. The 
objective transmission amplitude can be accurately mod- 
eled by a Gaussian function: 



T ob j(0) =T A exp 



f 2 sin 2 
4^2 



(6) 



where Ta is the axial transmittance amplitude, and £ 
is a length scale characterizing the objective for a given 
wavelength. This effect can readily be incorporated in 
terms of an effective reduced beam waist w e ^ given by 
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(7) 



We represent the optical force by an efficiency factor 
in the usual way [4]: 



x ^m,i(^i) Jm-i(kpsinQ) 



where d J 



,{0) are the matrix elements of finite rota- 
tions [22], and J m are the Bessel functions of integer or- 
der. The parameter 7 = f/w is the ratio of the focal 
length to the beam waist at the entrance aperture of the 
objective (not at the sample). It determines the fraction 
of available beam power that fills the objective aperture 
(filling factor [4|). The fraction of available power trans- 
mitted into the sample chamber is given by 



Q 



n\P/c 



(8) 



where P is the local laser power at the sample and c is the 
velocity of light. The evaluation of the trapping force in 
terms of a partial- wave series is similar to that described 
in [151 ] . We derive two separate contributions: Q = Q s + 
Q e . Q e represents the rate of removal of momentum from 
the incident beam. Its axial component is 
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2 ' ds s exp(— 27 2 s 2 ) 



y/(l ~ S 2 ) (N 2 - S 2 ) 



Vn 2 



- s z ) 
(5) 



where j ranges from one to infinity and m ranges from 
— j to j, the star denotes complex conjugation, aj and bj 
are the Mie coefficients [24J , Gj ?m is given by Eq. ((4]) and 
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dO sin Vcos cos 1 exp(-7 2 sin 2 0) T{6) d j m 1 (9 1 ) J m -i(kp sin 0) . 



(10) 



The transverse components in cylindrical coordinates are 
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G± m = J dO sin0vW#sin#i cxp(- 7 2 sin 2 6) T{0) d? m .±i.i( e i) J m _i(A^sin (9) e' 



(13) 



Q s represents minus the rate of momentum transfer 
to the scattered field. Its components are given in the 
Appendix. 

In this paper, we are particularly interested in the ax- 
ial component (in order to derive the trapping threshold 
and the stable equilibrium positions, if any), and in the 
transverse trapping stiffness 



n\P dQ p 
c dp 



where the derivative is taken at the equilibrium position 
z eq , which lies along the z-axis (p = 0). To obtain k±, we 
first perform termwise differentiation of the partial- wave 
series for Q p and then perform the summation numeri- 
cally. Similarly, z eq is computed numerically as the root 
of Q z (z eq ) = 0. When multiple equilibrium positions are 
found, we take the most stable root, corresponding to the 
deepest potential well. 

The distance L between the interface and the parax- 
ial focus in Eq. (|3]) is not experimentally accessible. In 
the experiments, the (inverted) microscope objective is 
first moved down until the trapped microsphere just 
touches the interface, and then the desired height in 
the sample chamber is obtained by moving up the ob- 
jective through a known distance d (for additional ex- 
perimental details, see next section). In order to mimic 
the experimental conditions, we adopt the following pro- 
cedure. We first compute the critical distance L c for 
which the equilibrium position is such that the sphere 
touches the interface, by numerically solving the equa- 
tion Q z (z eq — a — L c ) = for L c . This determines the 
paraxial focal plane for the initial configuration. By mov- 
ing up the objective through a distance the paraxial 
focal plane is displaced by Nd. Accordingly, we evalu- 
ate the equilibrium position and the transverse stiffness 
taking 



unit local power as functions of microsphere radius and 
height in the sample chamber. Aiming at absolute cali- 
bration, we propose to make a blind comparison between 
theory and experiment, with no fitting parameter. Hence 
the relevant trap parameters, such as the power P at the 
sample and the beam waist w at the entrance aperture 
of the objective, were measured directly, and the result- 
ing values were plugged into the MDSA model. When- 
ever possible, two different techniques were employed for 
measuring each parameter, and the results were checked 
against each other for consistency. The experiments were 
performed independently at UFMG lab and COPEA lab, 
reproducibly at intervals of several months. 

The two labs employed very different setups, represen- 
tative of those most often found in practice. Fig. [2] is 
a schematic drawing of the Diode Laser Setup (DLS), 
employed at UFMG lab. 
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III. EXPERIMENTAL PROCEDURES 



FIG. 2: 
(DLS). 



Schematic representation of the diode laser setup 



A. Experimental setups 

The most important force calibration data for practical 
applications are the values of the transverse stiffness per 



A Nikon TE300 inverted optical microscope (infinity 
corrected) with 100X N.A. 1.4 Plan Apo CFI objec- 
tive is employed for trapping, microsphere observation 
and scattered light collection. To one port is attached 
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a CCD camera (CCD-72 DAGE-MTI) for visualization; 
in the other port, we use a photodetector (EGG - Pho- 
ton Counting Module, SPCM-200-PQ-F500), with col- 
lection diameter of 150 /im, mounted in Newport XY 
stages for precise positioning. The photodetector deliv- 
ers TTL pulses ready to be fed into a Brookhaven BI- 
9000AT digital correlator. An infrared (IR) diode laser 
(SDL, 5422-H1) operating at 832nm is used for trapping. 
A He-Ne laser (SP-127), operating at 632. 8nm, is the 
scattering probe. A 20nm width line filter is placed in 
front of the photodetector to eliminate IR and any light 
other than that of the He-Ne laser. A half- wave plate 
and polarizers are used to control the intensity and po- 
larization of the He-Ne incident and scattered light. A 
motor (m) is connected to mirror Ml, which drives the 
IR beam onto the objective. The purpose of this mo- 
tor is to move the IR beam and, consequently, move the 
trapped bead in relation to the fixed He-Ne laser beam, 
to obtain the backscattering profile. By determining the 
backscattering profile and measuring the backscattered 
light intensity autocorrelation function (ACF), one can 
obtain the decay time of the Brownian position fluctua- 
tions, and finally get the trap stiffness. This procedure 
is described in detail by Viana et al. [251 ]. 

The YAG Laser Setup (YLS), employed at COPE A 
lab, is shown in Fig.[3l The YAG laser beam is expanded 




CCD 



FIG. 3: Schematic representation of the YAG laser setup 
(YLS). 

by the lenses Li and L2 to a diameter of 10mm. The mi- 
croscope and objective are the same model as in the DLS 
setup. The expanded beam is led to the objective by mir- 
rors Mi and M2. Before the microscope epi-illumination 
port a diaphragm D is employed to control the beam 
radius. The microscope stage is moved by step motors 
in the x and y directions, and a Prior motion controller 



is used to get precise positioning (1 /im accuracy). We 
use a CCD camera (Hamamatsu C2400) to observe the 
samples. In front of the CCD is placed a line filter that 
blocks light at the wavelength 1.064 /im, to prevent scat- 
tered and reflected laser light from entering the view field. 
The signal from the CCD is fed to a Scion Digital Frame 
Grabber. The digitized images are analyzed with the 
NIH ImageJ program. 



B. Measurement of beam waist 

To test MDSA theory, it is not necessary to measure 
the waist at the sample region, which is not an input 
parameter. On the other hand, the beam waist at the 
entrance aperture of the objective (or more generally the 
beam transverse intensity profile) is an important theo- 
retical ingredient, which was measured independently by 
two different techniques: the CCD and the diaphragm 
methods. 

In the CCD method, the objective of the microscope 
is removed and the sample is replaced by a mirror. After 
taking off the line filter, the reflected beam is imaged by 
the CCD. To get a ruler calibration, we first increase the 
laser power so as to saturate the image and obtain a good 
definition of the diaphragm border, and then we mea- 
sure the corresponding diaphragm aperture (see Fig. [3]). 
Fig. [H shows an image of the expanded YAG beam (10 
mm diameter) and Fig. [5] shows the radial intensity pro- 
file obtained from Fig. [4j The intensity profile is fitted 
to a Gaussian function, 

I(p) = I e-&. (14) 

The beam waist radius is w — 2<j. 

From the fit we get <tyls = 2.3 ± 0.2 mm for the YAG 
beam. The error bar was obtained from the statistics 
of five experiments performed over a period of one year. 
The same method was applied to the diode laser beam, 
yielding <jdls = 1.1 ± 0.1 mm. 




FIG. 4: YAG laser beam image. 



In the diaphragm method, we vary the diaphragm 
aperture R and measure, as a function of i?, the transmit- 
ted laser power P(R). With the Gaussian beam profile 
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C. Measurement of power at the objective focus 
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FIG. 5: Beam waist measurement with the CCD method: 
YAG beam intensity variation in the transverse plane. The 
Gaussian fitting (solid line) yields ctyls = 2.3 ± 0.2mm. 



fit, we have 

P(R) = Pt (l-eSy (15) 

where P t is the total power of the beam incident on the 
diaphragm. 

Fig. [6] is a plot of P(R) for the YAG laser. Adjust- 
ing the measured values to (jTSJ) , we found a beam half- 
waist ctyls — 2.1 =b 0.2 mm. Here again, the error bar 
was obtained from the statistics of five experiments per- 
formed over a period of one year. Applying the same 
method to the diode laser we found a beam half waist 
&dls — 1.3 ± 0.1mm. Therefore, the two methods agree 
to within the error bars. 



To measure local power at the sample for the DLS 
setup, two different methods were used: the dual objec- 
tive method [26[ and a mercury microbolometer method 

Fig. [3 schematizes the dual objective method. On top 
of the microscope stage, with the inverted objective be- 
low, is mounted a second, identical objective, in the up- 
right position. Pe is the power at the entrance of the 
inverted objective, and P ou t is the power transmitted by 
the compound system. The transmittance is assumed the 
same for both objectives. The objectives are positioned 
using three Newport actuators and the microscope stage 
in order to get a collimated beam emerging from the sec- 
ond objective, coaxial with the beam entering the first 
objective. 

p , 

out 



Objective 2 



7^==^=? on 




Objective 1 



P E 



FIG. 7: Schematic representation of the dual objective 
method. 




R(mm) 

FIG. 6: Beam waist measurement with the diaphragm 
method: variation of transmitted power versus diaphragm ra- 
dius R for the YAG beam. 



From Pe and P out we obtain a value for the overall 
transmittance of the compound system. However, care 
must be taken, since the objective transmittance in the 
infrared is usually not uniform: it is smaller for rays at 
larger distances from the axisfcorresponding to larger 
angles in the sample region) [23[. Thus, the root-mean- 
square overall transmittance overestimates the average 
objective transmittance and the power at the sample (see 
Ref. [23] for details). 

The Gaussian transmission amplitude model (|6|) yields 
the transmittance function (with p = f sin 6) 

Wp) = tAexp (-7^2) > ( 16 ) 

where t& = T\ represents the axial transmittance [28[. 

Table [I] shows the measured values for the parame- 
ters defining the Gaussian transmittance function ([T6]h 
for the YLS and DLS laser wavelengths. The entrance 
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TABLE I: Parameters characterizing the Gaussian transmit- 
tance function given by Eq. JB}, for two different wavelengths. 
tA is the axial transmittance, and £ is the transverse length 
scale associated to the transmittance radial variation. 









f (mm) 


DLS (A 


0.832 /mi) 


0.53 ±0.01 


2.2 ±0.1 


YLS (A 


1.064 /xm) 


0.31 ±0.01 


3.3 ±0.1 



aperture radii of the objectives were 3.5 mm for both 
the DLS and the YLS. From these values and the laser 
beam waists at the entrance port, we computed the av- 
erage objective transmittance t for each setup. We found 
t = 0.20 ± 0.02 for the YLS, and t = 0.45 ± 0.03 for the 
DLS. As expected, both values are below the respective 
axial transmittances shown in Table H To evaluate the 
error bars, we considered the propagation of the errors 
in the measured beam half waists. 

The value for the DLS setup was checked with the 
help of a mercury microbolometer method [27[. The mi- 
crobolometer, inserted into the DLS setup, consists of a 
standard Corning microscope glass slide (thickness 170 
/im), with an O-ring of 1cm diameter and 0.5cm height 
glued onto it, filled with water and containing mercury 
droplets with sizes in the /im range. Since we use an oil 
immersion objective, the microbolometer does not have 
any glass-air interface. Moreover, the oil refractive in- 
dex (1.496) being very close to the glass refractive index 
(1.51), we can ignore refraction at this interface in the 
analysis. Since the typical radii of the mercury droplets 
are in the micrometer range, each droplet may be mod- 
elled as a sphere embedded within two semi-infinite me- 
dia (Fig. Ej. 




FIG. 8: Schematic representation of the microbolometer. 

The incident laser beam is focused by the objective 
onto the mercury bead, heating it. A steady-state situa- 
tion is achieved in a time scale of the order of 1 second. 
The mercury heat conductivity is about 13 times larger 
than the water conductivity. Therefore, temperature at 
the surface of the mercury droplet remains constant for 
a given laser power. The laser power is then slowly in- 
creased until the mercury bead jumps. This occurs when 
its surface temperature reaches the boiling point of water 
(T5 = 97.2°C at the laboratory), making the determina- 
tion of this point very easy and accurate. Knowing the 



mercury droplet radius, the heat conductivities of the 
medium (^i = 0.6791 W/mK for 100°C) and glass slide 
(^2 = 9.43 W/mK for 100°C), the absorption coefficient 
of mercury for the laser wavelength used (Abs = 0.272), 
one can determine the local laser power P. The objec- 
tive transmittance is obtained by measuring the mercury 
droplet radius as a function of the input power producing 
the 'jump' (see for details). We found t = 0.40±0.04, 
in agreement with the dual-objective result within error 
bars. 



D. Measurement of bead radius 

A thorough test of MDSA theory requires covering a 
broad range of microsphere radii. Precision calibrated 
commercial beads are available only for a restricted size 
range. The DLS setup employs uncalibrated immersion 
oil droplets in de-ionized water, produced by a sonicator, 
yielding an almost continuous range of sizes, from sub- 
micrometers to tens of micrometers. The radius of the 
droplets employed for testing needs to be directly mea- 
sured. 

Up to a few /im, diffraction precludes reliable size de- 
termination by videomicroscopy. To measure radii in the 
range from 0.5 to 2.5 /im, we employed an alternative 
method, based on the properties of the free Brownian 
motion performed by the beads when the optical trap is 
turned off. We infer the hydrodynamic drag coefficient f3 
from the diffusion coefficient D by using the Einstein re- 
lation, and then obtain the radius from Faxen's extension 
of Stokes's law [2[. 

To check the method, we first applied it to a deionized 
water solution of polystyrene spheres of radius a = 1.52 ± 
0.03/im. A 10 _4 % solution is placed within a chamber 
built from an O-ring with ~ 1 cm diameter and ~ 0.3 cm 
width, glued on a coverslip with wax-candle, and covered 
by another coverslip to avoid evaporation. 

To control the bead height h = z ± L (see Fig. [1]), 
we first trap it with the optical tweezers and, employing 
the microscope knob, we move it towards the coverslip. 
When the bead touches the coverslip we see a change in 
its image. This is our reference height (h = a) and we use 
it to place the bead at a desired height by acting again 
on the microscope knob. This procedure allows height 
determination with an uncertainty of 0.5 /im. 

After positioning the bead, we choose a small image 
area (~ 50 pixels x ~ 50 pixels), in order to get a good 
frame capture rate, 27 to 28 frames per second. Then the 
position p = p p((j)) of the bead center-of-mass on the xy 
plane (the three-dimensional position is r = p ± zz) is 
measured at every 1/28 s. We turn the optical potential 
periodically on (~ 0.5 s) and off (~ 0.5 s) by shutting the 
beam. In Fig. [9] we show a typical result for the radial 
distance p(t) as a function of time. Note that p = 0, 
corresponding to the location of the trapping beam, also 
represents the trapping position on the xy plane. 
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FIG. 9: Radial distance to the beam axis as a function of 
time when the optical trap is periodically switched on and 
off. When the trap is off, the bead undergoes free Brownian 
motion. 



The mean square displacement is given by 

{{5p)\t))=ADt, (17) 

where the diffusion coefficient is given, for temperature T, 
by the Einstein relation (ks = Boltzmann constant) [30] 



D 



k B T 

(3 



(18) 



The drag coefficient f3 (for motion parallel to the glass 
interface) is given by Faxen's extension of the Stokes law 
in terms of the bead radius a and the height h [29[: 



6irr]a 



JL 

16 
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45 
256 



\h) 16 U/ 



(19) 



From the data illustrated in Fig. [9j we calculate ((dp) 2 ), 
excluding the regions where Sp < 0.01/im, the precision 
limit for determination of the bead center of mass using 
usual centroid-finding algorithms 0. Such regions corre- 
spond to the time intervals with the trap turned on. The 
time interval t in Eq. ([T7]) is an integer multiple of the 
inverse frame acquisition rate St. 

In view of the time translation symmetry of free Brow- 
nian motion, we may combine data from different time 
intervals, allowing us to improve our statistics in comput- 
ing ((Sp) 2 ). The data in Fig. [9] represents an ensemble 
of which each time window with the trap switched off is 
a realization. For each realization, we compute the aver- 
age mean square displacements ((Sp) 2 ) between nearest 
neighbors (separated by St) , between next nearest neigh- 
bors (separated by 2St), and so on, according to the for- 
mula 



N—k 



((6p) 2 (kSt)) 



N 



{P[U + k)5t] - P [j5t}} 2 , (20) 



where N is the total number of frames in a given time 
window and p is the radial position vector of the bead 
center of mass on the xy plane (note that only p = \p\ is 
plotted in Fig. [9]). 

The final result (((Sp) 2 (kSt))), the ensemble average 
over all time windows, is plotted in Fig. [10] as a function 
of time t = kSt. From the linear fitting to ([17]) . we obtain 
the diffusion coefficient D. 
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t(s) 

FIG. 10: Mean square bead displacement as a function of 
time. 

From the experimental value for D, we evaluate the 
drag coefficient (3 using Eq. JT5J). This procedure is re- 
peated at various heights. Fig. [11] shows the resulting 
values of /3 as a function of h. 
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FIG. 11: Measured drag coefficient f3 as a function of bead 
height h. The nominal bead radius is (1.52±0.02) /im, whereas 
the value measured from the Brownian motion is (1.5=b0.2) /xm 

The continuous line in Fig. [TT] represents the curve fit 
to ([T9|) . Near the coverslip wall, a small variation in h 
generates a substantial variation [cf. ([T9]) )] in the corre- 
sponding value of /3. To account for this source of error 
in height determination, we introduce an off-set ho in the 
curve fitting function. From the fit we get ho = 0.3pm 
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and a = 1.5 ± 0.2/im. The latter is in good agreement 
with the nominal radius 1.52 ± 0.02/im of the calibrated 
bead. 

The dashed lines in Fig. [11] correspond to the calcu- 
lated values of (3 for a variation of ±10% in the sphere 
radius a. We see that the data points are spread out 
between the two dashed lines, showing that the sphere 
radius is measured with 10% uncertainty. For droplets 
with a > 2.5/im, the radius was measured by videomi- 
croscopy and ruler calibration. Fig. [12] shows the image 
of a polystyrene sphere of measured radius 5.8 ± 0.1/ira 
and Fig. [13] shows the corresponding grey-level profile. 
We take the sphere diameter to be the distance between 
the centers of the sigmoid branches in Fig. [13l and the 
error bar as 1 pixel (~ 0.1 /im). 




FIG. 12: Bead image: radius equal to 5.8/xm, scale bar 10/im. 
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FIG. 13: Radial grey level of the bead image shown in Fig. 
1121 The bead diameter, defined by the distance between the 
side hills, equals 11.6 /im. 

For the YLS setup, we employed a few available cali- 
brated beads in the range below 2 /im, and a polydisper- 
sion of uncalibrated beads above this value, with radii 
measured by videomicroscopy. 



E. Measurement of trap stiffness 

For the DLS, trap stiffness is measured by analyzing 
Brownian motion in the optical potential well (as op- 
posed to the free Brownian motion employed for bead 
size measurement) (25|. A He-Ne laser is used to probe 
the bead motion, by detecting the intensity fluctuations 
of backscattered light. By measuring their time intensity 
autocorrelation function (ACF) with a digital correlator, 
one obtains the Brownian relaxation time of the bead 
(see Fig. Ell . 




0.98 I I 

0.05 0.1 0.15 0.2 0.25 

t(s) 

FIG. 14: Typical normalized ACF fitted to Eq. (|2l) . Inset: 
same ACF in a semilog plot. 

After normalization, the ACF can be fitted to the fol- 
lowing equation: 

g(t) = 1 + A± exp(-t/r±) + A z exp(-t/r z ), (21) 

where A± and A z are the amplitudes and tj_ and r z are 
the Brownian motion decay times in the radial and ax- 
ial directions, respectively. From the resulting value of 
t_l, we get the transverse stiffness from (see Ref. [25| for 
details) 

«L = — , (22) 

where the drag coefficient (3 is known from the measure- 
ment of the bead radii (see Sec. III.D). 

The trap stiffness in the YLS was measured using 
Faxen's law and videomicroscopy. After trapping a bead 
and positioning it at the desired height, the microscope 
stage is set to move laterally with a known velocity v. The 
Stokes force displaces the bead to a new off-axis equilib- 
rium position. The entire process is recorded, digitizing 
the images with a Scion frame grabber. The procedure 
is repeated until one has 5 to 10 different values for the 
stage velocity in each direction. Velocities are taken small 
enough to probe only the harmonic sector of the potential 
well (linear regime). The images are analyzed in order 
to infer the displacement from the equilibrium position 
5p eq as a function of v, checking that the dependence is 
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linear, of the form Sp eq = av. The trap stiffness is then 
obtained from the angular coefficient a by the relation 

«L = (23) 
a 



IV. RESULTS 

Following the procedures discussed in Sec. Ill, we have 
measured the ratio k±/P as a function of a, at different 
heights h above the glass interface. Measurements of all 
relevant data allowed for a direct comparison with MDSA 
theory, with no adjustable parameters. In addition to the 
measurements of beam waist and bead radii discussed 
in Sec. Ill, we also measured the refractive indices of 
immersion oil (712 = 1.496) and deionized water (ni = 
1.343) at the diode laser wavelength Ao = 832 nm by the 
minimum deviation method. 

For the polystyrene beads employed in the YLS, we 
have taken the value ri2 = 1.576 reported in Ref. [3l| 
for the YAG laser wavelength Ao = 1064 nm, and we 
measured the refractive index of water to be n\ = 1.332 
at this wavelength. Much less important is the refractive 
index of the glass slide, which is relevant for the spherical 
aberration effect, but not for Mie scattering: n = 1.51. 

For the DLS, the objective vertical displacement d was 
adjusted to have (approximately) the same height h for 
all bead sizes. If the equilibrium position of the bead 
with respect to the paraxial focus were independent of L 
(distance between focus and glass wall), we would have 
a height 

h = a + Nd (24) 

after displacing the objective by a distance d. A different 
displacement d was taken for each a, so as to keep a con- 
stant h, as given by ([24]) . Because the optical potential 
depends on the L, the actual height differs from ([M]) by 
an amount < a. This difference is negligible with regard 
to the Faxen correction of Stokes law (see Section III.D), 
but it might be important for a quantitative evaluation 
of the optical force and stiffness. Thus, it is taken into 
account when computing the force and the trap stiffness. 

However, the dependence with L and h turns out to 
be negligible for the DLS, because the laser beam under- 
fills the objective entrance aperture. Thus, the objective 
produces a more paraxial beam in the sample chamber, 
rendering the aberration effect irrelevant. Indeed, the 
laser beam waist is rather small (odls = 1.1 =t 0.1mm, 
see Sec. III.B) in this setup. Moreover, a diaphragm 
(radius r c = 2.0 mm) was employed, cutting off the tails 
of the transverse intensity profile. This translates into a 
reduced effective numerical aperture NA e ff = nsin# e ff, 
where sin^eff = r c /f, where / is the objective focal 
length. In terms of the entrance aperture radius r e , the 
sine condition also yields / = r e / sin#o? where #0 rep- 
resents the semi-aperture angle for plane-wave illumina- 



tion. Therefore, 

NA eff = - NA, (25) 

r e 

where NA = 1.4 is the objective numerical aperture. For 
our objective, we have r e = 3.5 mm, yielding NA e ff = 0.8, 
which corresponds to # e ff = 32°. Thus, we replace #0 by 
this value in (|4|) and ([5]) when computing the theoretical 
values for the DLS. 

For the YLS, on the other hand, the aberration ef- 
fect is important, since the laser beam overfills the en- 
trance aperture. Plane- wave components at angles > 
arcsin(A^) are refracted into evanescent waves in the sam- 
ple chamber. Since we perform the experiments at a dis- 
tance of several wavelengths from the glass interface, the 
contribution of evanescent waves is neglected. When us- 
ing the MDSA model presented in Sec. II, we replace 
#0 by arcsin(A/") = 61.9°, corresponding to an effective 
numerical aperture NA e ff = m = 1.332. 

As discussed in Sec. II, the radial variation of the ob- 
jective transmittance is taken into account by using the 
effective waist defined by (0). From the values mea- 
sured for the beam waist (Sec. III.B) and the length 
£ characterizing the transmittance variation (Table I of 
Sec. III.C), we find w e R = 2.2 mm and iu e ff = 3.5 mm for 
the DLS and the YLS, respectively. 



A. Trap stiffness 

Fig. [15] displays the results for the DLS at 
h = (3.1 ± 0.5) /im. Data points are averages of 4 inde- 
pendent measurements (each one lasting 100s) for ev- 
ery microsphere; vertical error bars show the associated 
standard deviations. Full and dashed lines correspond to 
MDSA and GO (geometrical optic, i. e., WKB) theo- 
ries, respectively. For 0.04 /im < a < 0.52 /im, the former 
predicts that no stable trapping is possible, because no 
equilibrium position is found in this range [32]. This is 
indicated by the vertical dotted line break. Correspond- 
ingly, experimental points cluster around different values 
in the neighborhood of the threshold at a = 0.52 /im. 
In this range, the microspheres often escaped from the 
trap during the measurement interval. Scattered data 
points closely below threshold arise from microspheres 
that stayed in the trap for at least 3 measurements 
(around 300 s). Since the harmonic potential well ap- 
proximation breaks down in this metastable region, such 
data should be regarded mainly as confirmation of un- 
stable below-threshold behavior. 

As shown analytically in Ref. [15[, the GO values cor- 
respond to size averages over the oscillations of the exact 
curve for sufficiently large values of the size parameter 
27rnia/Ao. Fig. [T5l shows that GO still provides not too 
bad estimates of stiffness for sizes down to a ~ Ao- This is 
not surprising, since incident rays that contribute when 
the sphere center is axially located are unaffected by the 
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FIG. 15: Transverse trap stiffness (divided by the local power) 
as a function of bead radius for oil beads (712 = 1.496) in 
deionized water (m = 1.343) at height h = 3.1 /im, for the 
DLS (Ao = 0.832 /xm, objective underfilling). Circles: exper- 
imental points (with error bars). Solid line: MDSA theory. 
Dashed line: geometrical optics. The threshold for trapping 
is indicated by a dotted vertical line. 



dominant diffraction and resonance effects in Mie scat- 
tering, which are most important for rays incident near 
the edge of the microsphere, and WKB is a good approx- 
imation down to size parameters of order unity outside 
of this edge domain H. Prior estimates of the domain of 
applicability of GO [2] placed it beyond size parameter 
30 (a > 5Ao ). Our present theoretical and experimental 
extension of that domain has obvious practical relevance. 

Fig. [16] shows the results for the YLS setup with 
d = (3.0 ± 0.5) /im, which corresponds to sphere heights 
h — a^ (2.7 ±0.5) /im according to Eq. ([M]) . Polydisper- 
sions allowed taking many data points for larger a, but 
only a few (calibrated) polystyrene sizes were available 
for smaller a. For such microspheres, at least 4 indepen- 
dent measurements were taken (vertical bars represent 
the corresponding standard deviations). Here we find 
(both experimentally and theoretically) stable trapping 
below the peak of the stiffness curve. As compared to 
the DLS, the major difference is the use of (moderate) 
objective overfilling in this setup. 

The amplified scale in the inset of Fig. [16] reveals the 
persistence of oscillations around the GO curve in the 
tail of the theoretical curve. They arise from interference 
among multiple reflections at the upper and lower mi- 
crosphere interfaces. As shown in [ll|, these oscillations 
can be derived analytically from the partial- wave series. 
The corresponding period is given by Aa = \Xo/n2 ~ 
0.169 /im, in very good agreement with the solid line 
curve in the inset of Fig. [16] For the DLS setup, on 
the other hand, the oscillations are distorted (Fig. H5J) 
by the beam aperture constraint. 

Our measurement accuracy is not sufficient to ver- 
ify the presence of the interference oscillations, a very 
demanding experimental challenge. However, the inset 
shows that theory already predicts a band of halfwidth 
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FIG. 16: Same conventions as in Fig. [15] for polystyrene 
beads (712 = 1.576) in deionized water (m = 1.332) at 
d — 3 /im, for the YLS (Ao = 1.064 /xm, moderate objective 
overfilling). Inset: amplification of the region 2.8 /im < a < 
4.0 /im. 



of the order of 5%, about half our measurement uncer- 
tainty, around the GO average, within which data points 
are expected to fall. 

For the lowest measured a, theory predicts two equi- 
librium points, but only the most stable one, closer to 
the paraxial focus, is considered. We have checked the 
numerical results against the Rayleigh limit (very small 
a). In this regime, the force is proportional to the in- 
tensity gradient of the incident field alone, allowing us 
to compare our results with those of [33[. Spherical 
aberration leads to the appearance of two local inten- 
sity maxima along the axis, which correspond to stable 
equilibrium positions in the Rayleigh limit. Outside the 
Rayleigh regime, but with small values of a, there are 
still multiple equilibria. This is clearly connected with 
the spherical aberration effect (thus explaining why no 
such effect takes place in the DLS). 

We may enhance the effect of spherical aberration by 
increasing the vertical displacement d of the objective. In 
Fig. [T71 we show the YLS results for d = 15 /im, corre- 
sponding to h — a « (13.2 ±0.5) /im. The inset shows that 
the period of oscillation for large a is not modified by the 
aberration effect, as expected. In the Rayleigh regime, on 
the other hand, we find seven stable equilibrium points. 
They correspond to local intensity maxima lying along 
the axis, between the glass slide and the paraxial fo- 
cus, that result from interference fringes bordering the 
focal region in the diffraction theory of spherical aber- 
ration [33[. As we decrease a, starting from the value 
a = 0.44 /im, the number of equilibrium points rapidly 
increases and then saturates in the Rayleigh regime, as 
shown in Fig. [18] In this regime of multiple equilibria, it 
is very difficult to achieve trapping and to measure the 
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stiffness. Moreover, the potential well around the equi- 
librium point near the paraxial focus is usually very shal- 
low in this case (see for instance the solid line in Fig. [19] 
for a — 0.265 /im). Thus, spherical aberration degrades 
trapping efficiency for small a, producing a 'threshold' for 
stable trapping at a = 0.44 /im. This is indicated by the 
vertical dotted line in Fig. [TT] (further discussion of mul- 
tiple equilibrium points is given in Sec. IV. B). By com- 
paring this figure with Fig.[16l we conclude that spherical 
aberration decreases the peak stiffness value and slightly 
displaces the peak towards larger values of a. On the 
other hand, the effect on stiffness is very small for larger 
a values (GO limit). 

Increasing the height above the glass slide does not 
change the stiffness in the DLS, as expected for a more 
paraxial beam. Measured and computed stiffness results 
at h = (8.6 ± 0.5) /im (not shown) differ very little from 
those of Fig. [151 
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FIG. 17: Same as Fig. [H at d = 15 /im. 
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FIG. 18: Number of stable equilibrium points as a function 
of sphere radius, for the YLS with d — 15 /im. 



B. Bead hops 

As discussed in Sec. IV. A, the optical potential for 
small sphere radii contains several equilibrium points be- 
cause of the spherical aberration effect. In the course 
of Brownian motion of the bead, it may hop from some 
local minimum of the optical potential to a more stable 
neighboring minimum. We have observed this effect and 
measured the size of the hop. 

In Fig. [191 we pl°t the optical potential along the z-axis 
in units of /c#T as a function of z/a, for three different 
values of the objective vertical displacement: d = 15 /im 
(solid line), d = 6.5 /im (dashed line) and d = 3.0 /im 
(dotted line). We take parameters corresponding to the 
YLS, with a 0.265 /im and P = 20 mW. 

Besides the equilibrium point nearest to the parax- 
ial focus (z = 0), Fig. ITUI shows additional equilibrium 
points, located between z = and the glass slide. They 
decrease in number and in degree of stability (as mea- 
sured by the corresponding well depth) as d decreases. 

The arrow in Fig. [T9l indicates the most stable equilib- 
rium position when d = 15 /im. As d decreases to 6.5 /im, 
its well depth remains close to 60/c#T, and then starts 
to decrease much faster below this point (this local min- 
imum also approaches the paraxial focus, as expected, 
because the interval between this point and the glass 
slide, which contains all local intensity maxima, becomes 
shorter as d decreases). At d = 6.5 /im, the potential well 
near the focus is already much deeper, and the distance 
between the two equilibrium points is 2.1 /im. 




FIG. 19: Optical potential along the z-axis (in units of UbT) 
for the YLS, with a = 0.265 /im and local power P = 20mW. 
Objective displacement: d = 15 /im (solid line), d = 6.5 /im 
(dashed line) and d = 3.0 /im (dotted line). The arrow points 
to the most stable equilibrium point for d = 15 /im. 

These predictions were tested by video microscopy, fol- 
lowing the image of the trapped bead as the objective 
was displaced downwards (thus bringing the paraxial fo- 
cus closer to the glass slide), from d = 15 /im down to 
d = 3.0 jam. The image sharpness depends essentially on 
the bead position with respect to the paraxial focus. Ini- 
tially, the bead is trapped at the position indicated by 
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the arrow in Fig. [19l As we bring down the objective, at 
first the image does not change, indicating that the bead 
follows the displacement of the paraxial focus, keeping an 
approximately constant distance to this point ('tweezer' 
effect). Hence, in this first stage the bead remains in 
the same potential well. However, at d = 6.5 /im, we ob- 
served a sudden change in the image (see Fig. |20|) . Since 
the bead size is a fraction of the wavelength, the image 
behaves approximately like the diffraction pattern of a 
point source by the circular objective aperture. On the 
focal plane (bottom image), this approaches the standard 
Airy pattern. Away from the focal plane (top image), it 
corresponds to a Fresnel diffraction pattern, with central 
spot brightness (in this case, dark) related to the number 
of Fresnel zones within the aperture [34| . 

In order to calibrate the bead-focus distance, we re- 
peated the experiment with the bead attached to the 
glass slide. In this case, the distance between bead and 
paraxial focus is known for each position of the objective. 
By comparing the resulting images with those obtained 
with the optically trapped bead, we derive the length of 
the (upward) hop to be 2.2 ± 0.5/im, in good agreement 
with the MDSA prediction for the distance between the 
two equilibrium points at d = 6.5 /im, shown in Fig. [19] 
(dashed line). 
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FIG. 20: Bead hop in the optical potential, a) Schematic 
representation of the position of the trapped bead, paraxial 
focus F (origin of the z axis) and glass slide S for d > 6.5/xra 
(top, bead below F) and d < 6.5/ira (bottom, bead close to 
F). b) Bead images for d > 6.5/xra (top) and d < 6.5/ira 
(bottom), c) Calibration of bead-focus distance with bead 
attached to the glass slide. 

When the objective is moved upwards, we do not ob- 
serve a reverse hop. Instead, the initially trapped bead 
is lost as d approaches 15 /im. This is also consistent 
with the solid line in Fig. [T9j as the well near the focus 
becomes shallower, the bead cannot climb the optical po- 
tential barrier and reach the deeper well below the focus 
(indicated by the arrow in Fig. [T9|). This is in line with 
our previous discussion about the difficulty of trapping 
small beads when the paraxial focus is too far from the 
glass slide ('threshold' of Fig. [T7j). 



V. CONCLUSIONS 

We have demonstrated that the MDSA ab initio the- 
ory of trapping forces in optical tweezers, formulated 
solely in terms of experimentally accessible parameters, 
already comes close to achieving a basis for absolute cal- 
ibration. Indeed, we have verified in detail its conse- 
quences concerning trap stiffness and trapping thresholds 
for two very different setups, peak locations, relationship 
with and domain of validity of geometrical optics, size 
and height dependence resulting from interface spherical 
aberration, covering the whole range from the Rayleigh 
regime to the geometrical optic one, and including the 
description of multiple equilibria and 'hop' effects. 

Among many attempts (Sec. I), it is the only treat- 
ment that satisfies requirement (iv) of Sec. I, asymptot- 
ically approaching geometrical optics in the mean. This 
is specially important in view of our finding that geo- 
metrical optics is already a reasonable approximation, as 
far as trapping stiffness is concerned, for radii of the or- 
der of the wavelength (a consequence of the absence of 
near-edge diffraction effects). 

The combined experimental procedures and experi- 
mental results reported here, with well-defined and ac- 
curately measured parameters, far outnumber all previ- 
ously available optical tweezers stiffness data. Compar- 
ison with MDSA theory shows that, within the typical 
10% order of magnitude of the error bars, theory and 
experiment are generally in very good agreement, for a 
broad size range and two very different (and often em- 
ployed) setups, one with underfilling and the other one 
with moderate overfilling of the objective. The theory 
accurately predicts the trapping threshold and the loca- 
tion of the stiffness peak, which is extremely sensitive 
to beam shape, as well as the effects of spherical aber- 
ration arising from refraction at the interface, including 
multiple equilibria situations. 

The largest deviation from theoretical predictions oc- 
curs within the stiffness peak, where the steeply chang- 
ing curve is most sensitive to additional perturbations, 
for the YAG laser setup. As the microsphere radius de- 
creases in this region, one enters the domain of multi- 
ple equilibria and Brownian fluctuations are enhanced. 
Note (Fig. [T6|) that the experimental result lies below 
the MDSA prediction, indicating that missing effects pre- 
sumably arise from degradation of the focal region. Nat- 
ural candidates are objective aberrations in the infrared 
[effect (i) in Sec. I]. In fact, agreement is better for the 
diode laser setup, which employs a more paraxial beam 
and a wavelength closer to the visible range. 

Other effects not taken into account in MDSA theory 
include reverberation (multiple beam reflections between 
the bead and the glass slide) and the effects of evanes- 
cent waves beyond the critical angle, as well as possible 
bead surface distortions or contaminations. In applica- 
tions requiring accurate force measurements, in view of 
these perturbations, it seems advisable to stay at least a 
couple of wavelengths away from the glass slide. Account- 
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ing for the deviations may yield new insights concerning 
bead-surface interactions at close range. 

Polarization effects were not discussed in the present 
work: we found the transverse stiffness to be practically 
independent of polarization. Polarization effects might 
be relevant, on the other hand, when the sphere is at 
a distance from the axis larger than or of the order of 
the radius [35[. They will be discussed in a forthcoming 
work. 

It has sometimes been argued that absolute calibra- 
tion is an impossible aim, in view of the multiplicity of 
effects that need to be taken into account, as well as the 
errors in parameter determination. We do not share this 
pessimistic outlook. Mie scattering, the basic interaction 
involved, is understood and verified at a level of precision 
approaching that of quantum electrodynamics |3[ . While 
we had to employ uncalibrated microspheres to cover a 
wide size range, NIST-traceable calibrated microspheres 
are commercially available. Achieving an understanding 
and control over additional perturbing effects can turn 
out to be relevant not only to improving the performance 
of conventional optical tweezers, but may have potential 
applications to a variety of new techniques in precision 
microscopy. 
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